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We study the five-orbital Hubbard-Holstein model for iron pnictides with small electron-phonon 
interaction due to Fe-ion Einstein oscillators. Using the fluctuation-exchange (FLEX) approxi- 
mation, orbital fluctuations evolve inversely proportional to the temperature, and therefore the 
resistivity shows linear or convex T-dependence for wide range of temperatures. We also analyze 
the Eliashberg gap equation, and show that s-wave superconducting state without sign reversal 
(s+_l_-wave state) emerges when the orbital fluctuations dominate the spin fluctuations. When both 
fluctuations are comparable, their competition gives rise to a nodal s-wave state. The present study 
offers us a unified explanation for both the normal and superconducting states. 
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The many-body electronic states and the pairing mech- 
anism in iron pnictides have been significant open prob- 
lems. By taking account of the Coulomb interaction and 
the nesting of the Fermi surfaces (FSs) in Fig.l (a), fully- 
gapped sign-reversing s-wave state (s±-wave state) had 
been proposed [1, 2]. Experimentally, both T c and an- 
tifcrro (AF) spin correlation increases as x decreases in 
BaFe2(Asi_a;P K )2 [3]. In contrast, T c in LaFcAsOi-^F.,; 
at x = 0.14 increases from 23 K to 43 K by applying 
the pressure, whereas AF spin correlation is almost un- 
changed [4]. Thus, the relationship between T c and spin 
fluctuation strength seems to depend on compounds. 

On the other hand, orbital-fluctuation-mediated s- 
wave state without sign reversal (s ++ -wave state) 
had been proposed based on the five-orbital Hubbard- 
Holstein (HH) model [5, 6], which is the Hubbard model 
introduced in Ref. [1] with the addition of the electron- 
phonon (e-ph) interaction term due to Fe-ion Einstein 
oscillations. Within the random-phase-approximation 
(RPA), it was found that rf-orbital fluctuation is induced 
by small e-ph interaction, not by the Coulomb interac- 
tion alone. Especially, empirical relation between T c and 
the As-Fe-As bond angle (Lee plot) [7] is naturally ex- 
plained. Recently, theoretically predicated orbital fluc- 
tuations [5, 6] had been detected via the substantial soft- 
ening of shear modulus [8]. The s ++ -wave state is consis- 
tent with the robustness of T c against randomness [9-11] 
as well as the "resonance-like" peak structure in the neu- 
tron inelastic scattering [12]. 

However, spin/orbital fluctuations obtained by the 
RPA are reduced by the self-energy correction. There- 
fore, to confirm the orbital fluctuation scenario, it is de- 
sired to analyze the many-body electronic states beyond 
the RPA. For this purpose, the fluctuation-exchange 
(FLEX) approximation [13] would be appropriate in that 
the absence of spin/orbital order in 2D systems, known 
as Mermin- Wagner theorem, is rigorously satisfied [14]. 

In this letter, we analyze the five-orbital HH model 



for iron pnictides using the FLEX approximation. In 
the normal state, large orbital fluctuations induce highly 
anisotropic quasiparticle lifetime on the FSs as well as the 
T- linear or T-convex resistivity p [15-17]. The large or- 
bital fluctuations also introduce the s+ +-wave supercon- 
ducting (SC) state for wide range of parameters, and the 
competition between orbital and spin fluctuations lead to 
the nodal s-wave state. We propose that the orbital fluc- 
tuation is the origin of both the s++-wave SC state and 
the non-Fermi-liquid-like behavior in the normal state. 

In the FLEX approximation [13], the 5x5 self-energy 
matrix £ in the orbital representation is given by 

= ^EE^w 3l4 M G ^( fc -9). (i) 

q l 2 h 

where N is the number of k meshes, and we denote 
k = (fe, e n ) with fermion Matsubara frequency e n = 
(2n + l)7rT, and q = (g, uj n ) with uj n = 2mrT. Here, 
U represents the Z 2 , XZ, YZ, X 2 - Y 2 and XY or- 
bitals, which are hereafter denoted as 1, 2, 3, 4 and 5, 
respectively [1]. G is the 5x5 Green function matrix in 
the orbital basis, and V 12 is the 25 x 25 interaction term 
for the self-energy given as [18] 

- I ( f* _ f c )x irr (g)(f s - f c ) + r]r s + \t c (2) 

where the irreducible susceptibility is given by 

til^iM = -^J2 G hh(k + q)G hh (k), (3) 

k 

and the spin (orbital) susceptibility is obtained as x s ^ — 
X s ™ (1 - f s ( c )x irr ) _1 - Here, f s = S (f c = -C - 2V(uj n )) 
is the irreducible vertex for spin (charge) channel; S and 
C represent the Coulomb interaction in the multiorbital 
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model introduced in Refs. [1, 5, 6, 18]; Their matrix ele- 
ments consist of the intra-orbital Coulomb U, the inter- 
orbital Coulomb U', Hund's coupling J and the pair hop- 
ping J'. Here we assume that J = J' and U = U' + 2J. 

The vertex V(w n ) in T c represents the electron-electron 
(el-el) interaction mediated by e-ph interaction. For ex- 
ample, we show non-zero Vj;/ ;mm / for I, V , m, ml = 2, 3, 4 
in Fig. 1 (b), where g(ui n ) = gw^/{uil + cj£) is pro- 
portional to the phonon Green function; g — g(0) is the 
effective el-el interaction for u> n — 0, and ojd is the De- 
bye frequency [5]. Other than Fig. 1 (b), V has many 
non-zero off-diagonal elements as explained in Ref. [6], 
since the Fe-ion oscillation (non-A lg mode) induces var- 
ious inter-orbital transitions. This fact gives rise to the 
prominent orbital fluctuations at low frequencies, while 
the charge susceptibility x c (<?) = 12i m Xu, mm (q) is not 
enhanced due to the cancellation [5, 6]. In the present 
study, we drop ladder-type diagrams by g(u) n ), which is 
justified when wd <C Ef [5, 6]. For the same reason, 
g{uj n ) is absent in V s . 

In the FLEX approximation, we obtain G and S self- 
consistently using the Dyson equation G _1 = (G°) _1 — S. 
In multiband systems, the FSs are modified from the orig- 
inal FSs due to the self-energy correction. To escape from 
this difficulty, we subtract the constant term [S(fe, +i0) + 
£(fe, — i0)]/2 from the original self-energy, correspond- 
ing to the elimination of double-counting terms between 
LDA and FLEX [19]. Hereafter, we fix J/U = 1/6, 
cjd = 0.02eV, and the electron filling n — 6.1 except 
for Fig. 4. Because of the smallness of the FSs in Fig. 1, 
fine k meshes are required for a quantitative study. We 
take N = 128 x 128 fc meshes that is four times of that 
used in Ref. [19], and 1024 Matsubara frequencies. Then, 
we obtain reliable numerical results for T > O.OleV. 
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FIG. 1: (color online) (a) FSs in the unfolded zone. The dot- 
ted circles represent the cold-spot given by the orbital fluc- 
tuation theory. The cold-spot is composed of xz/yz-oihitals. 
(b) Phonon-mediated el-el interaction (V) for 2, 3, 4 orbitals. 

We begin with the electronic property in the normal 
state. Hereafter, the unit of energy is eV. First, we dis- 
cuss the quasiparticle damping rate 7& on each FS, which 
is given by the imaginary part of the self-energy in the 
band-diagonal representation. In Fig. 2 (a), 7^ repre- 
sents the damping due to spin (orbital) fluctuations for 
T = 0.015 and U — 1.8, which is given by substituting 
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FIG. 2: (color online) (a) fc-dependence of induced by 
the spin (orbital) fluctuations on each FS. (b) T-dependence 
of S s ( c ) = (1 — a s ( c )) _1 . (c),(d) T-dependence of p. p = 
1 corresponds to (ha c )/e 2 ~ SOO/icmfl when the interlayer 
distance is a c — 0.6nm. 



V s = |f s x s r s (if c x c f c ) in eq. (1). The horizonal 
axis is the azimuth angle for k point with the origin at 
r(M) point for FS1,2 (FS4). The relation 7fc 7^ + 7^ 
is satisfied since the third term in eq. (2) is very small. 
We will see below that the value U = 1.8 can reproduce 
moderate AF spin fluctuations observed in e-doped com- 
pounds, and it is consistent with U ~ 2 reported by x-ray 
absorption spectroscopy (XAS) [20]. 

In the case of g(0) = g = 0, the relation 7 s 3> 7 C holds 
since the orbital fluctuation is very small. However, 7 C 
gradually develops with g, and 7 s ~ 7 C at g — 0.26. In 
Fig. 2 (a), 7 C 3> 7 s at g — 0.3; the corresponding dimen- 
sionless coupling is just A = gN(0) ~ 0.2. In contrast, 
7 s decreases with g, due to the suppression of x s by 7°. 

As shown in Fig. 2 (a), the momentum dependence of 
7! on each FS is small, although AF spin correlation is 
well developed. In contrast, 7^ on FS4 (e-pockets) takes 
the minimum value at 8 ~ for g = 0.3; this "cold- 
spot" is important for the transport phenomena. Since 
the cold spot is on the e-pocket, the Hall coefficient Rn 
will be negative, which is consistent with experiments 
[15, 16, 21]. In the case of high- T c cuprates, various non- 
Fermi- liquid-like transport phenomena (e.g., violation of 
Kohler's rule) originate from the cold/hot spot structure 
as well as the backflow (=current vertex correction) due 
to spin fluctuations [22]. Therefore, appearance of the 
cold spot in Fig. 2 (a) indicates that the orbital fluc- 
tuations are the origin of striking non-Fermi-liquid-like 
transport phenomena in iron pnictides [15, 16, 21]. 

In Fig. 2 (b), we show how the orbital and spin 
fluctuations develop as T decreases: In the FLEX, the 
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spin (orbital) susceptibility is enhanced by the spin (or- 
bital) Stoner enhancement factor S s / C -) = (1 — a s ( c )) , 
where a s r c -\ is the maximum of the largest eigenvalue of 

f s ( c ^x 1TT (q,0) with respect to q. Then, a 8tC = 1 corre- 
sponds to the spin/orbital order, although it is prohibited 
in 2D systems by the Mermin- Wagner theorem [14]. In 
case of U = 1.8 and g — 0, large S s (> 10) is induced 
at q w Q = (tt,0) (i.e., x s (QiO) x S s )- S s gradually 
increases as T drops, which is a typical critical behavior 
near the AF magnetic quantum-critical-point (QCP) [23]. 
When g > 0, x c (<7> 0) is enhanced at q — and q = Q 
almost equivalent ly [6]. At g = 0.3, large S c (3> 10) 
is induced at both q « Q and 0, and it increases ap- 
proximately proportional to T" 1 . Thus, it is confirmed 
that both ferro- and AF-orbital fluctuations show critical 
evolusions near the orbital QCP. 

Next, we discuss the resistivity p due to the orbital 
and spin fluctuations. By neglecting the backflow, the 
conductivity is obtained by 

1 ^ f°° dta ( df(u)\ | x _ , . n .,a , A , 

k,a J ~ co V 7 

where a is the band index, f(co) is the Fermi distribution 
function, v x a k is the velocity of band a, and Gk_ a (uj + i0) 
is the retarded Green function for band a in the FLEX 
approximation. Figure 2 (c) and (d) show the obtained 
the resistivity p — l/a xx for U = 1.2 and 1.8: In case of 
U = 1.2, p shows a conventional sublinear (concave) In- 
dependence at g = 0. p increases with g due to the orbital 
fluctuations, and almost T-linear resistivity is realized at 
g = 0.2. At g = 0.25, p shows a superlinear (convex) 
T-dependence. In case of U = 1.8, p is linear-in-T at 
g = 0, while it shows a clear superlinear T-dependence 
at g = 0.3. Thus, we stress that non-Fermi-liquid resis- 
tivity is not a direct evidence for the spin fluctuations. 
In inFeAsO compounds, T c increases as the radius of 
ranthanide ion Ln 3+ decreases, and the T-dependence of 
p changes from concave to convex [17]. This experimen- 
tal correlation between T c and p(T) may be understood 
in terms of the orbital fluctuation scenario. 

Next, we discuss the SC state. In the presence of dilute 
impurities (n; mp <C 1), the linearized Eliashberg equation 
in the orbital basis is [5]: 

X E A w (k) =- ^ E W ^,mAk-k')G' mim2 {k') 

k' : rrii 

x A m2m3 (k>)G' mim3 (-k>) + <^,(e„),(5) 

where Aw (fc) is the gap function and Ae is the eigenvalue 
that reaches unity at T = T c . SSigma represents the 
impurity-induced gap function. (G") _1 = (G) _1 — 5T, n , 
where G is the Green function given by eq. (1), and <J£" 
is the impurity-induced normal self-energy. The pairing 
interaction W in eq. (5) is 

w(q) = lr s x s (q)r s Jr c x c (<z)r c + If* - If c , (6) 
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FIG. 3: (color online) (a) U-g phase diagram given by solving 
the linearized Eliashberg equation at T = 0.015. Nodal s- 
wave gap state is obtained in the shaded area for ni mp = 0.02 
and 0.05, and solid lines (dotted lines) represent the boundary 
between fully-gapped s+_-wave (s++-wave) state. Dashed- 
dotted line denotes a c — 0.98. (b) s++-wave gap (A_e = 0.59) 
for U = and g = 0.24. (c) s+_-wave gap (A B = 0.49) for 
U — 1.8 and g = and. (d) Nodal s-wave gap (A_e = 0.28) 
for U = 1.2 and g = 0.15. 

where x s ' c is given by the FLEX approximation for 
"-imp = 0, because of the fact that the fully self-consistent 
FLEX with impurity-induced self-energy leads to un- 
physical reduction in \ s > unless vertex correction is taken 
into account [22]. The first (second) term in eq. (6) works 
to set A FS i i2 • A FS3: 4 < (> 0). 

In the T-matrix approximation, (5S n,a is given as 

tf£y(e„) = ni mp Tij(e n ), (7) 
tf£?-(e n ) = n imp ^Tii(e n )fi m (e n )T jm (-e n ), (8) 

lm 

where Ty(e„) = 1(1 — 7g(e„)) _1 is the T-matrix in 

the normal state [9]; g(e n ) = ^J2k^ k ( €n ) 1S ^ ne 
local normal Green function, and I is the local im- 
purity potential that is diagonal in the orbital basis. 
We put 1 = 1 hereafter. In eq. (8), fij(e n ) = 
7fJ2k,im G n( k )^-irn(k)G jm (-k) is the linearized local 
anomalous Green function. 

In Fig. 3 (a), we show the U-g phase diagram ob- 
tained by the FLEX approximation. The dashed-dotted 
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line represents the condition a c = 0.98 at T = 0.015, 
corresponding to g — 0.25 ~ 0.3. (In the RPA, the same 
condition is satisfied for g = 0.21 ~ 0.23 [6].) Therefore, 
prominent orbital fluctuations emerge for A = gN(0) < 
0.2 even if the self-energy correction is taken into ac- 
count. On the other hand, a s = 0.95 (0.92) for U = 1.8 
and g = (0.3) in the FLEX approximation, although 
Uqi = 1-25 for a s = 1 in the RPA. Thus, the renormal- 
ization in a s is rather larger than that in a c , because of 
the difference in the coefficients (in factor 3) between the 
first and the second terms in eq. (2). 

Next, we solve eq. (5) with high accuracy using the 
Lanczos method at T — 0.015. Then, s ++ -wave gap 
function is obtained around the line a c = 0.98; Figure 
3 (b) shows the s++-wave gap for g — 0.24 and U = 
(Xe = 0.59). On the other hand, s±-wave gap is ob- 
tained when g is sufficiently small; Figure 3 (c) shows 
the s±-wave gap for U — 1.8 and g — (Xe = 0.49). 
When riimp = 0, the gap function changes from (b) to 
(c) discontinuously on the phase boundary in Fig. 3 (a), 
as found in Ref. [6]. When n imp > 0.02, however, gap 
function changes smoothly during the crossover. Then, 
line-nodes inevitably appear on FS3,4 in the shaded area 
in Fig. 3 (a) ; Figure 3 (d) shows the nodal s-wave gap for 
U = 1.2, g = 0.15 and n imp = 0.02 (A B = 0.28). Thus, 
both regions for s ++ -wave and nodal s-wave states are 
extended by small amount of impurities, although Xe for 
the latter state is reduced by impurities. A nodal s-wave 
solution at rii mp = with larger Xe may be obtained by 
considering a three-dimensional (3D) nodal-line structure 
in 3D tight-binding model [24]. 

Here, we discuss that line nodes originate from the 
competition between the orbital and spin fluctuations: 
The electrons at 9 ~ {tt/2) on FS4 is composed of 
orbital 2,3 (4). Since the orbital 4 is absent in FS1,2, 
the nesting-driven spin correlation between the orbital 
2,3 on FS1,2 and the orbital 4 on FS3,4 is weak. (That 
is, X24 42(Q) *^ X22 22( < 3)-) On the other hand, both 
X24 42( < z) an d X22 22(?) are developed here [5, 6]. There- 
fore, when orbital and spin fluctuations are comparable, 
Apsi,2 • Aps4 is negative (positive) at 9 ~ (7r/2) due to 
the orbital-dependence on the spin correlation. 

In Fig. 4 (a) , we show the filling dependence of Ae for 
the s++-wave state (U — 1.2, g = 0.25), and that for the 
s+_-wave state (U = 1.8, g = 0). We note that FS1,2 
disappear for n > 6.3. The value of Ae for the s++-wave 
state decreases monotonically with n, while Ae for the 

s_| -wave state is rather insensitive to n, maybe because 

the temperature, T — 0.02, is rather high. Figure 4 (b) 
shows that Ae for the s++-wave state (U — 1.2, n = 6.1) 
increases with g. 

In summary, we performed the FLEX approximation in 
the multiorbital HH model for iron pnictides. It was con- 
firmed that orbital-fluctuation-mediatcd s++-wave state 
is realized by small e-ph interaction g. As increasing 
the value of g, both the T c of s ++ -wave state and the 
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FIG. 4: (color online) (a) n dependence of Ae for s++- and 
s+_-wave states at T = 0.02 and ni mp = 0. (b) g dependence 
of Ae for s ++ -wave state at T = 0.015 and n lmp — 0. 



resistivity p are increased, and the latter changes from 
T-concave to T-convex. This correlation between T c and 
p is consistent with experiment [17]. The region of s++- 
wave or nodal s-wave states is enlarged in the presence 
of small amount of impurities. Thus, the present or- 
bital fluctuation scenario presents a unified explanation 
for both normal and SC electronic states. 
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